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fi I We discuss the non-linear evolution of the angular momentum L acquired by pro- 

■ tostructures, like protogalaxies and protoclusters, due to tidal interactions with the 

surrounding matter inhomogencitics. The primordial density distribution is assumed 
to be Gaussian and the non-linear dynamics of the collisionlcss mass fluid is followed 
using Lagrangian perturbation theory. For a Cold Dark Matter spectrum, the inclu- 
sion of the leading-order Lagrangian correction terms results in a value of the rms 
ensemble average (L 2 ) 1 / 2 which is only a factor of 1.3 higher than the corresponding 
linear estimate, irrespective of the scale. Consequently, the predictions of linear theory 

\ are rather accurate in quantifying the evolution of the angular momentum of proto- 

■ structures before collapse sets in. In the Einstein-de Sitter universe, the initial torque 
' is a good estimate for the tidal torque over the whole period during which the object 

is spun up. 
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1 INTRODUCTION 



The problem of the acquisition of angular momentum by protostructures in the universe is of considerable interest in theories 
. ' of galaxy and cluster formation. A widely accepted view is that the present luminous structures acquired their spin via 
■ gravitational tidal interactions with the surrounding matter inhomogeneities (Hoyle 1949; Peebles 1969; Doroshkevich 1970; 
i-j 1 White 1984; Barnes & Efstathiou 1987; Hoffman 1986, 1988; Heavens & Peacock 1988; Ryden 1988; Quinn & Binney 1992; 
Eiscnstein & Loeb 1995; Catelan & Theuns 1996). So far, however, the theoretical analysis of the growth of the tidal galaxy 
angular momentum has been essentially limited to the linear regime, during which the galaxy spin grows proportionally to 
the cosmic time t (Doroshkevich 1970; White 1984). 

In this paper we examine analytically, for the first time, the question of how the galaxy tidal angular momentum evolves 
during the mildly non-linear regime. Previous attempts in this direction may be found in Peebles (1969) and White (1984). 
However, in contrast to their approach, we employ actual perturbative solutions of the dynamical equations that describe the 
motion of the fluid. More in detail, we apply Lagrangian perturbation theory. 

The Lagrangian approach shows to be ideal in treating the evolution of the galaxy spin, because it is powerful in describing 
the non-linear growth of the mass-density fluctuations on one hand (Zel'dovich 1970a, b; Buchert 1992; Bouchet et al. 1992; 
Catelan 1995); and on the other hand the usual difficulty of inverting the mapping from Lagrangian coordinates q to Eulerian 
coordinates x is completely by-passed. This is because the angular momentum L is invariant with respect to the Eulerian 
or Lagrangian description. 

The layout of this paper is as follows: in the next section we first briefly review the basics of the Lagrangian theory 
and the perturbative solutions of the Lagrangian fluid equations. Next, we compute within this framework the perturbative 
corrections to the linear tidal angular momentum iJ- 1 ' acquired by a protoobject. The resulting expressions are then simplified 
by calculating their averages over the ensemble of realisations of the linear Gaussian gravitational potential ip^ 1 ' for objects 
with given inertia tensor. We then compare the non-linear spin corrections with the results of the linear analysis recently 
performed by Catelan & Theuns (1996). In the main text we restrict ourselves mainly to the case of a flat universe, leaving 
the more involved treatment of closed and open universes to Appendices. 
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2 NON-LINEAR EVOLUTION OF THE TIDAL ANGULAR MOMENTUM 

We consider a Newtonian pressureless and irrotational self-gravitating fluid embedded in an expanding universe with arbitrary 
density parameter £1 but zero cosmological constant. Such a fluid is assumed to mimic the behaviour of matter on scales smaller 
than the horizon. Furthermore, we assume that luminous objects like galaxies and clusters of galaxies started to grow due to 
gravitational instability around primordial positive density fluctuations S in this fluid. 

We indicate by x comoving Eulerian coordinates, from which physical distances may be obtained according to the law 
r = a(t)x, where a(t) is the expansion scale factor and t the standard cosmic time. We will use the temporal coordinate r of 
which the differential is defined by 

dr = aT 2 dt , (1) 

instead of t, since this allows a considerable simplification of the formalism when dealing with the gravitational interactions 
in a generic non-flat Friedmann universe (Shandarin 1980). The peculiar velocity and the peculiar acceleration simplify when 
using r instead of t to: 

dx 

— = x = o(r) u(x, r) , (2) 

=X=g(x,T) . ( 3 ) 

The dimensionless time r is negative, as discussed in Shandarin (1980), and the initial cosmological singularity at t = 
corresponds to r = — oo. In the open models, the infinity of the cosmic time, t = +oo, corresponds to r = — 1; in the Einstein- 
de Sitter universe, t — +oo corresponds to r = 0; and in the closed models the contraction phase starts at r = 0. In terms of 
the density parameter Q, one has 

t = -\/^F(l-fi)- 1/2 , (4) 

where k is the curvature constant (k = — 1 for open universes and k = 1 for closed universes). The case Q — 1 (k = 0) is a 
singular point for the transformation (4) and in this case we take r = — (3t) -1 ^ 3 , which corresponds to using a(t) = (3t) 2//3 
or to = t/a 3 / 2 = 1/3, which defines the unit of time. The scale factor a(r) may then be written for all Friedmann models as 
o(t) = (r 2 + k)-\ 

The linear evolution of angular momentum of protoobjects is most easily analysed using the Zel'dovich (1970a, b) 
formulation (see White 1984; Catelan & Theuns 1996) and the mildly non-linear spin growth is most easily analysed using the 
Lagrangian perturbation theory. We recall that the Zel'dovich approximation coincides with the linear Lagrangian description. 

We will essentially adopt the formulation of the Lagrangian gravitational theory for a collisionless Newtonian fluid 
presented in Catelan (1995; see also references therein) but note that in the present paper, the variable r has opposite sign 
and the growth factor of the linear density perturbation is normalised differently. An alternative formulation of the Lagrangian 
theory may be found in Buchert (1992). 



2.1 Basic tools: Lagrangian theory 

In Lagrangian formulation, the departure of the mass elements from the initial position q is described in terms of the 
displacement vector S, 

x(q,r)sq+S(q,r). (5) 

The trajectory x(q, t) of the fluid element originally at q satisfies the Lagrangian 'irrotationality' condition and the Poisson 
equation given by (Catelan 1995) 

(-a/3 1 X ya = , (6) 

x%'x 0a = a{r){J - 1] , (7) 

respectively, where e Q /3 7 is the totally antisymmetric Levi-Civita tensor of rank three, £123 = 1, and summation over repeated 
Greek indices (where a = 1,2,3) is understood. In these equations, q(t) = 6a(r) and J = 1/(1 + <5) is the determinant of 
the Jacobian of the mapping x — > q. The determinant J is non zero until the first occurrence of shell-crossing (see, e.g., 
Shandarin & Zel'dovich 1989). Furthermore, x a p = dx a /dqp and denotes the cofactor of x a p: we recall that the latter 
is a quadratic function of x a p and consequently, the master equations (^j) and ([?]) are cubic in x a p. In general, x a p is not a 
symmetric tensor: x a p = xg a if, and only if, the Lagrangian motion is longitudinal in which case x can be obtained from the 
gradient of a potential. 

The irrotationality condition and the Lagrangian Poisson equation can be written in terms of the displacement field S as: 
£ a p~t [(1 + V ■ S) 8p a - Spa + S^ CT ] SV = , (8) 
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[(1 + V ■ S) S afJ - S afJ + S%] S 0a = a(r)[J(q, r) - 1] , (9) 

where V = V q and the symbol 5 a/ g indicates the Kronecker tensor. These equations are the closed set of general dynamical 
equations for the displacement vector S describing the motion of a collisionless fluid in the Lagrangian {q}- space, embedded 
in an arbitrary Friedmann universe and subject to the Newtonian gravitational interaction of the mass fluctuations J -1 — 1. 
We briefly summarise their perturbative solutions in the next subsection. 



2.2 Lagrangian perturbative solutions 

The exact Lagrangian equations (8) and (9) are non-linear and non-local in the displacement S (see the discussion in Kofman & 
Pogosyan f 995) and it is undoubtedly very difficult to solve them rigorously. A possible alternative is to seek for approximate 
solutions by expanding the trajectory S in a perturbative series, the leading term being the linear displacement which 
corresponds to the Zel'dovich approximation: S = Si + S 2 + S 3 + . . ., where S n = 0(S") is the n-th order approximation. 
Note that a perturbation series of this form needs to include at least the third-order term S3 to capture the essential physics 
contained in the cubic equations (8) and (9). 

For the sake of simplicity we will limit ourselves to the Einstein-de Sitter universe in the main text. In this universe, 
a — t~ 2 . The reader interested in more general Friedmann models is addressed to Appendix A. We neglect decreasing modes. 



2.2.1 First-order approximation: Zel'dovich approximation 

The first-order solution to equations (8) and (9) is separable in space and time and corresponds to the Zel'dovich approximation 
(Zel'dovich 1970a, b): 

Si(q,T) = D(T)S")(q)Ej)(r)V^»(q) ; (10) 

the function D(t) is the growth factor of linear density perturbations and is given in the Einstein-de Sitter universe by: 

D(t) = t~ 2 = a(r) . (11) 

The expression of D(t) valid in a generic Friedmann model is reported in Appendix A. The function t/>' 1 - ) (q) is the (initial) 
gravitational potential. For later use, we define its Fourier transform, i/^ 1 ' (p) = J dqip^(q) e _ipq , where p is the comoving 
Lagrangian wave vector. The Fourier transform of the linear density field, S^'(p,t) = D(r)Si(p), is related to V>^(p) via the 
Poisson equation, (p) = p~ 2 5i(p) . 



2.2.2 Second-order approximation 

The second-order solution is again separable with respect to the spatial and temporal variables and describes a longitudinal 
motion in Lagrangian space: 

S 2 (q, r) = E(t) S< 2 > (q) = E(r) W> (2) (q) . (f 2) 

The growing mode E(t) is for the Einstein-de Sitter model given by: 

E(r) = - V 4 . (13) 
The analytic expression for E(t) in a non-flat universe and the expression for -i/>' 2 '(p) are reported in Appendix A. 

2.2.3 Third-order approximation 

The third-order solution S3 corresponds to three separable modes, two longitudinal and one transverse, denoted by subscripts 
a, b and c, respectively: 

S 3 (q,r) = F„MSi 3, (q) + ft(r)S< ;i) (q)+F c (r)T< 3 )(q) = F a (r)V^(q) + F t MV*« 3, (q) + F c (r)Vx A< 3 >(q) . (14) 
The growing modes F a , F b and F c for a flat universe are, respectively, 

Mr) = -ir- 6 , (15) 
Fb{r) = +^r- 6 , (16) 
Fc(r) = -It- 6 . (17) 
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We stress that the transverse mode T ( } does not describe any physical vorticity in the fluid since the latter is assumed to be 
irrotational. Rather, the occurrence of a transverse component is due to the fact that a Lagrangian frame of reference is not 
inertial. Consequently, this 'fictitious' term is required to obtain a correct physical description of the motion. Unfortunately, 
this transverse mode is forgotten in some of the relevant literature on the subject (e.g. Bouchet et al. 1995). 

The general expressions for F a , Ft and F c for a non-flat universe as well as the potentials V>i 3 ' (p) and y^ 3 ' (p) and the 
transverse components T a are given in Appendix A. 

Comparing equations (11), (13) and (15-17), it is clear that the perturbative expansion for S is in fact a Taylor series in 
the variable D(t) — t~ 2 . As the general expressions reported in the Appendix A testify, this is no longer true in a non-flat 
universe. However, it is shown there that the higher-order growth factors can be approximated exceedingly well by powers of 
D: E oc D 2 and F oc D 3 , so the expansion is still 'close' to a Taylor expansion. We will consider next how the perturbative 
series in S translates into a perturbative series for the angular momentum L. 



2.3 Non-linear spin dynamics 

Let us consider some volume V of the Eulerian x-space. The Cartesian coordinate system is assumed to be centred at the 
centre of mass of V. Since we are interested in the intrinsic angular momentum of the mass contained in V, we disregard the 
centre of mass motion. 

The angular momentum L of the matter contained at the time t in the volume V is 

L(t) = p b (t) a(t) 4 I dx[l + «5(x, t)] x(t) X u(x, t) = Vo [ dq [q + S(q, r)] X dS ^' T ) , (18) 
Jv(t) Jr dr 

where we have substituted the time variable r in favour of the standard cosmic time t in the second integral. Here, the matter 
density field is p = pt(l + S), where pb(r) is the background mean density and S is the density fluctuation field, which is 
assumed to be initially Gaussian distributed and r)o = a 3 p b ; the peculiar velocity field is denoted by u (see, e.g., Peebles 1980). 

The second equality in equation (18) stresses the important fact that the integral over the Eulerian volume V may 
be written equally well as an integral over the corresponding (initial) Lagrangian volume F. This enables us to apply the 
Lagrangian description of the Newtonian gravity previously reviewed. The linear regime (Zel'dovich approximation) has 
been fully analysed in this way by Doroshkevich (1970), White (1984) and Catelan & Theuns (1996), whereas its Eulerian 
counterpart was studied by Heavens and Peacock (1988). We can extend the Lagrangian analysis of the evolution of the 
angular momentum L(r) to the non-linear regime by applying perturbation theory to equation (18). Perturbative corrections 
to S(q, r) (Bouchet et al. 1992; Buchert 1994; Catelan 1995 and references therein) then give perturbative corrections to L(r): 

h=0 

where we have defined 

L<">(T) = £ no [ dq S,(q,r) X dS »~^ , (20) 

with So = q, hence L' -* = 0. Since we will be interested in calculating second-order corrections to the ensemble average (L 2 ), 
we need to compute corrections to L up to third-order. After briefly reviewing the results of the linear theory, we summarise 
the final expressions of the corrections L^ 2 - 1 and L^K The reader interested in the details of the calculations is addressed to 
the Appendix B. 



L(r) 



= V °l 



dq 



£si(q,T) 



dr 



J2 s ^- r ) 



.3=0 



2.3.1 Linear approximation 

The linear Lagrangian theory corresponds to the Zel'dovich approximation and the first-order term in equation (19) is given 
by: 

L (1) (r) =nob(j) j dqqx W W (q) . (21) 

If ipW(q) is adequately represented in the volume F by the first three terms of the Taylor series about the origin q = 0, then 
each component L^\t) may be written in a compact form as (White 1984): 

La ] {t ) = D(t ) e al3j V ( pl J al , (22) 
where we introduced the deformation tensor at the origin: 
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V$=V${0)=d fl d ir 1,W{0), (23) 
and the inertia tensor of the mass contained in the volume F 

Ja-i = r)o J dq q a g 7 . (24) 

Equation (22) shows that the linear angular momentum L/ 1 -* is in general non-zero because the principal axes of the 
inertia tensor J a p, which depends only on the (irregular) shape of the volume F, are not aligned with the principal axes of 
the deformation tensor V^l , which depends on the location of neighbour matter fluctuations. The temporal growth of tidal 
angular momentum is completely contained in the function D(r), which behaves as D(r) — — 2t~ 3 ~ t in the Einstein-de 
Sitter universe, as first noted by Doroshkevich (1970). Finally, if F is a spherical Lagrangian volume, then ~ e Q( 3 7 = 0. 
Consequently, the matter contained initially in a spherical volume does not gain any tidal spin during the linear regime (see 
also the discussion in White 1984) . 

2.3.2 Second-order approximation 

The second-order term in equation (19) involves the second-order displacement S^ 2 - 1 : 

L (2) (r) = r !o J dq q X ^ = 770 £(r)|dqqx V^ (2) (q) . (25) 

Note that the second-order term which follows from the product Si x dSi/dr is identically zero since is involves the product 
Vi/> < - 1 - > x VV'^- The potential ip 1 ' 2 ' 1 is determined by the potential ip^ through the equation (64) in Appendix A. Note that, 
since E oc r~ 4 , one has E cx t~ 5 hence the second-order terms grows oc i 5//3 in the Einstein-de Sitter universe. This growth 
rate was first derived by Peebles (1969). If we represent t/>' 2 - ) (q) in F by the first three terms of a Taylor series around q = 0, 
as we did before for ip^, we obtain for the a-component: 

L^(r) = E(r)e a01 vfl J„ 7 , (26) 
where 

Vf^vf 1(0) = 8,8^(0). (27) 

We can call 2?^ the second-order deformation tensor. The component may be written in terms of the second-order shear 
tensor £ ' 2 2 as 

L% ) {T) = E{T)e olBl £f}j all (28) 
where 

4 2 j = V™ - \{V ■ S< 2 % CT = {d p d„ - i^ CT V 2 )V» (2) . (29) 

An analogous relation is valid for the linear term (see Catelan & Theuns 1996). Note that the non-linear dynamical evolution 
modifies only the deformation tensor and not the Lagrangian inertia tensor. Furthermore, if F is a sphere, then again 
This is in contrast to Eulerian perturbation theory, where the angular momentum of an Eulerian sphere does grow in second- 
order perturbation theory (Peebles 1969; White 1984). 

2.3.3 Third-order approximation 

The two longitudinal modes Si 3 ', and the transverse one T^ 3 ', and the coupling between first and second order displace- 
ments originate the following third-order spin corrections: 

L (3) (r) = m J dqq X ^ + m J dq (Sr X ^ + S 2 x ^) = L< 3 > (r) + L< 3) (r) + L< 3 > (r) + l/ 12 >(r) , (30) 
where 



H ' (r) = VO F h (r) J dq q X V^' (<0 . ( 31 ) 
L( 3) (r) = r ?0 F c (r)^dqqxT (3) (q) , (32) 



'r 

,7/3 



and h — a, b. Here, the temporal functions are such that F a oc Fb oc F c oc r oc t 1 in the Einstein-de Sitter universe 



Proceeding as in the previous cases by expanding the fields ipi 3 \ and in Taylor series around q = 0, one gets the 



expressions, 
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Lfl{r) = F c {r) e a p~, T$ Jp a . 

In these last equations, we defined the third-order deformation tensors 



h pa 



2?i 3 i(0)=a^ 3) (0) 



and again h — a, b and T£g = dT^ 1 /dqp. In addition, one has 



0), 



L (12) (r) = Vo [D(t)E(t) - E(t)D(t)] J dqVV (1> 



x Vi/> 



(2) 



where (DE — ED) oc r 7 oc t 7 ^ 3 in the Einstein-de Sitter universe. Taylor expanding as before results in 



Li 12) = (DE - ED) e af3j 



no r d fl ip m (o) d^> (o) + v% E# j ari 



(1) T)(2) 



A more explicit expression for the contribution L^ 12 ' in terms of integrals in Fourier space may be found in the Appendix B. 
Note that if one considers the volume F to be centered on a peak of the underlying density distribution, as is presumably 
appropriate when studying the formation of collapsing objects, the first term is zero since is an extremum at the origin 
in that case. 



(33) 
(34) 

(35) 
(36) 
(37) 



For a single collapsing region enclosed in a volume V it is enough to evaluate equation (19) at the time of maximum 
expansion tm to compute its angular momentum. After tm the angular momentum essentially stops growing since the collapsed 
object is less sensitive to tidal couplings (Peebles 1969). However it is more useful to compute the mean angular momentum 
of the object averaged over an ensemble of realisations of the gravitational potential random field tjj^ 1 ': this is particularly 
important in order to compare the theory against statistical results obtained from iV-body simulations or from observations. 
This programme is carried out in the next section. 



3 ENSEMBLE AVERAGES 

We simplify the previous results by considering the expectation value over the ensemble of realisations of ^i' 1 ' of the square of 
L, (L 2 )^ = (L 2 ), for objects with preselected inertia tensor. The underlying motivation for this is that it gives the appropriate 
analytical estimate to compare against numerical simulations that make use of the Hoffman-Ribak algorithm to set up a 
constrained density field that contains an object with given inertia tensor (Hoffman & Ribak 1991; van de Weygaert & 
Bertschinger 1996). Moreover, the resulting expectation value is still a good estimate for the angular momentum of Gaussian 
peaks in case where the correlation between gravitational potential field and inertia tensor - neglected here - can be properly 
taken into account, at least in the linear regime (see the discussion in Catelan & Theuns 1996): the exact calculation in the 
non-linear regime appears intractable analytically. In addition, the procedure could give some insight into perturbative spin 
corrections in the case of more generic primordial non-Gaussian statistics. 

Taking into account the mildly non-linear corrections, one has as leading terms: 

(L 2 ) = (L (1)2 ) + (L (2)2 ) + 2(L (1) ■ L (3) ) + o(r" n ) . (38) 

Note that the term (L' 1 ^ • l/ 2 '} is zero if the linear potential is assumed to be Gaussian distributed, as it involves an odd 
number of random fields. This term would represent the lowest -order perturbative correction to the linear term for the more 
general case of non-Gaussian statistics (work in progress). 



3.1 Perturbative corrections 

3.1.1 Linear approximation 

The linear ensemble average (L' 1 ' 2 } is computed and discussed extensively in Catelan & Theuns (1996): 

<L« 2 > = ^D(t) 2 (^ - 3 M2 ) o(R) 2 , (39) 

where the quantity o(R) 2 is the mass variance on the scale R, explicitly given by o(R) 2 = (27r 2 ) -1 dpp 6 P^(p)W(pR) 2 . In 
this expression, W denotes a filter applied to the input power spectrum to remove any ultra-violet divergence. We will use the 
Gaussian smoothing function W(pR) = exp(— p 2 R 2 /2). Equation (39) holds for any power spectrum (V>^(pi) (P2))y> = 
(27r) 3 <5d(pi + P2) Pii>{p) and the value of (l/ 1 ^ 2 ) depends on the normalisation of the spectrum; the symbol Sd indicates the 
Dirac function. The general expression (|39|) is independent of the details of the shape of the boundary surface of the volume 
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Y: it depends only on the quantities (ii and (12 , which are respectively the first and the second invariant of the inertia tensor 
J a 0- Specifically, denoting the eigenvalues of the inertia tensor by ti, 12 and 13, one has 

(ii = L1+L2 + (-3 (40) 
and 

M2 S Ht 2 + + t2t3 ■ (41) 

For a spherical volume, 11 = 12 = (.3, hence (jl\ — 3/i2 = 0, as we stressed before. For any volume Y one has (i\ — 3(i2 > 0. 
3.1.2 Higher-order approximation: (l/ 2 ^ 2 ) 

The calculation of the term (L^ 2 ) takes advantage of the results of the second-order approximation. The final expression, 
assuming a power-law spectrum Ps(p) = Ap" — P,pp 4 , may be written as: 

<L< 2 > 2 ) = A E(r f(p 2 - 3^2) E (2) (i?; n) , (42) 
where the function E^ 2 - 1 depends on the smoothing scale 7? and the normalisation of the spectrum A, as 

E< 2 > (R; n) = ^ J dpp n+2 [W(pR)] 2 A< 2 > (p; n) , (43) 



A<2) (Pi ! ") = dp2 P2+2 (l " O ' [(Pi P2 1 + P2 Pi 1 ? ~ 



40 2 l (n - 4)/2 . (44) 



Typically, these integrals have to be evaluated numerically. The results for more physical power spectrum, like the Cold Dark 
Matter (CDM) spectrum, are discussed in Appendix B, where the details of the derivation of the equation (42) are given as 
well. Surprisingly, the average (L/ 2 ' 2 ) factorises the same invariant (i 2 — 3(i2 of the inertia tensor J as appeared in the linear 
term (equation 39). 



3.1.3 Higher- order approximation: (i/ 1 ^ • L^, ) 

The calculation of the correction (I/ 1 )- L^ 3 '} takes advantage of the results of the third-order Lagrangian theory. The displace- 

(3) 

ment Sa does not induce on average any higher-order correction to the angular momentum, since it corresponds to radial 
motions of the fluid patches (see Appendix B for an explicit derivation), 

<L«.Li 3 >>=0. (45) 

Assuming again a scale-free power spectrum P^(p) = Ap n ~ 4 , the corrections due to the third-order displacements and 
T' 3 ^ are (see Appendix B for details): 



(L W. L («)) = ^D{T)F b {T)(n\-Z(i2)vf\R-n) , (46) 
2 



(L (1) -L< 3 >) = ~ D(t)F c (t)(» 2 - 3(12)^ (R;n) , (47) 
lo 



where the functions E^ and E^ are respectively 

Ef > (R; n) = -f^ / dpp n+3 [W(pR)} 2 A< 3) (p; n) , (48) 



E< 3 > (R; n) = -^7 / dpp n+A [W(pR)] 2 A< 3 > (p; n) , (49) 



(2-) 4 Jo 

and the integrands A^ and A^ 



f-OO P-\-l 

A.f\ P i-n) = J dp 2 p n 2 +1 j do (i~e 2 ) 2 [{pip- 2 1 +p 2 p- 1 1 )-29Y 1 , 

poo f~^~l 

Ai 3) ( P i;n)= J dp 2P n 2 J d6 (l-6 2 ) (6-P2PI 1 ) [(p lP 2 1 +P2Pi 1 )-2ey 1 . 

Again, these integrals have to be calculated numerically. The results for a more physical power spectrum are discussed in 
Appendix B. Note that once more these averages factorise out the invariant (i\ — 3(i2 of the inertia tensor. 



(51) 
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3.1.4 Higher-order approximation: (L"'> L/ 12 ') 

The perturbative correction (L^ ■ l/ 12 '} originates from the coupling between first- and second-order displacements. The 
calculation of this term is cumbersome and is detailed in Appendix B: we restrict ourselves to giving the final result: 

<L« L< 12 >) =^D(r) [D(r)E(r) - D(r)E(r)] (/i 2 — 3/i2 ) E* 12 ^ (R; n) , (52) 
where 

E (12) (i?;n) = ^-J ^^W(p 1 R)Wip2R)W(\p 1 -p2\R)P i ,(pi)P il (p2) fe^pffi Pi«(Pa-Pi)« [(Pi-PaJpL-pW*.] ,(53) 

where is the second-order kernel defined in Appendix A. This integral needs to be computed numerically both for scale- 
free and CDM spectra. Note that as previously this pertubative correction term to the spin is proportional to /if — 3/i2 . This 
remarkable property enables us to calculate the relative contribution of linear and higher-order spin terms. 

The time dependent growth factors of the various terms in equation (pq ) are illustrated in Fig. f and the results of 
numerically integrating the momentum contributions E for both power law and the CDM spectrum are shown in Fig. 2 as a 
function of the spectral index n. These wavevector integrations generally diverge at small wavelengths for scale-free spectra: to 
obtain finite expressions we have filtered out this ultra-violet divergence by smoothing the integrals over P2 (in equations 44, 
50 and 51) artificially with a Gaussian filter of width 0.5 x h~ x Mpc. The CDM integrals, in contrast, are finite. The numerical 
values of these momentum integrals will be used in the next section to estimate the relative contributions of the higher-order 
terms to the linear estimate of (L 2 ). The dependence of these relative momentum contributions on the smoothing radius R 
for a CDM spectrum is shown in Fig. 3, from which it is clear that, although the various E's depend strongly on R, the 
normalised contributions E/er 4 are more weakly dependent on the smoothing scale. 



4 ANGULAR MOMENTUM AT THE MAXIMUM EXPANSION TIME 

In this section we quantify the relative perturbative corrections to the linear angular momentum by computing: 



T (L 2 ) 1 _ ( L ' 2)2 ) , 9 (L (1) - L (l3) ) _ T (82) , T (13) , T (13) T (13) , r (112) , - , 

p 5 ) (LW 2 ) + a + " + c + 



where we recall that L/ 3 ' = L^ 3 ' +Lc 3 ' +L/ 12 - 1 hence L' 3 ' gives rise to four correction terms. From equation (45) we find 

immediately that T^ 13 ^ = 0. We compute the other correction terms at the time defined by D(tm)o~(M) — 1 on the mass scale 
M, which is close to the maximum expansion time found from extrapolating the spherical model (e.g., Peebles 1980). After 
tm, the protoobject starts collapsing and tidal torques are much less efficient in spinning up its matter content (Peebles 1969; 
Barnes & Efstathiou 1987). Therefore, we assume that the growth of the spin effectively ceases after maximum expansion of 
the object and identify the angular momentum at that time with the 'final' angular momentum. We understand that this is 
a partial description of the real world (see the discussion in Catelan & Theuns 1996). 
Collecting the expressions for the various corrections, we find: 

^ = (M ^pMrM)«(M)f = 0.22 (55) 



= -24^[ fl N"Wf = -o.i2 (57) 

D 2 D cr{Mp 

where the numerical values are calculated for the (flat) standard CDM model when filtered on the scale of 7? = 0.5 h~ x Mpc 
and at the maximum expansion time, i.e., when Da = 1. The factors E/DD, F b /D 2 D, F C /D 2 D and (DE - DE)/D 2 D do 
not depend on r for a flat universe; for a non-flat universe, their r dependence is extremely weak, in view of the excellent 
approximations E oc D 2 , Fj, iC oc D 3 and DE — DE oc D 3 , as illustrated in Fig. 5. Note that, with these approximations, the 
various T's do not depend on the normalisation of the power spectrum. 
We conclude that, since at maximum expansion 

(L 2 ) = (1 + T) <L (1)2 ) w 1.6 (L (1)2 ) , (59) 
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Figure 1. Time dependencies, for open, flat and closed universes, (top, middle and bottom curves respectively) (i/ 1 ) 2 ) 1 / 2 oc £>, 
(L (2) 2) l/2 K ^ ( L (l). L (3))l/2 (DF,) 172 , (L(l) • Li 3) >!/2 oc (_DF C ) 1/2 and <LM ■ l/ 12 )) 1 ^ k (d(DE - DE)) 1/2 . The latter 

four have been scaled by the indicated numerical factors to give the asymptotic behaviour (— r) — 5 for r — > — oo and are practically 
indistinguishable on the plot. The vertical line at r = — 1 denotes the infinity of physical time t for open universes. The term corresponding 
to DF a is not plotted since (I^ 1 ) ■ L^ 3) ) = 0. 

the linear estimate of (L 2 ) is roughly a factor 1.6 times smaller than the value obtained when taking into account the lowest- 
order non-linear corrections, for a standard CDM spectrum. Hence, */ (L 2 ) w 1.3 \/ (L/ 1 ) 2 }. From this we conclude that the 
predictions of linear theory are surprisingly accurate and that the dynamical perturbative corrections appear converged. 



5 SUMMARY AND CONCLUSIONS 

In this paper we analysed the growth of the tidal angular momentum L acquired by a protoobject (protogalaxy or protocluster) 
during the mildly non-linear evolution of the matter density perturbations, assuming the latter to be Gaussian distributed. 
The dynamics of the collisionless matter fluid is described using the Lagrangian approach in the formulation given by Catelan 
(1995). This formulation is very suitable to study the problem at hand, because the Lagrangian expressions are considerably 
simpler than their Eulerian counterparts, yet the protogalaxy's tidal spin is a vector invariant under the change of Eulerian 
to Lagrangian spatial coordinates, x and q respectively. Specifically, the difficult problem of inverting the mapping x = q + S, 
where S is the displacement vector, in order to recover the Eulerian quantities from the Lagrangian ones, is completely avoided. 

The strategy we follow is straightforward. The non-linear spin corrections \J- h \ where L' 1 - 1 is the linear angular momentum, 
are calculated approximating the fluid elements' trajectories S by the perturbative solutions S^ of the Lagrangian fluid 
equations (8) and (9). This leads to the expression (20). Since we are interested in computing the lowest-order perturbative 
corrections to the ensemble average (L* 1 - 12 ) for objects with given inertia tensor, we need to calculate corrections to L up to 
third-order. This has the added advantage that we take account of the full physical content of equations (8) and (9), since the 
latter are cubic in the displacement. The calculation is summarised as follows: from the knowledge of S = S1 + S2 + S3 (where Si 
corresponds to the displacement in Zel'dovich approximation), we deduce the corresponding corrections L = L 1 - 1 ^ +L' 2 - 1 -l-L^ 3 - 1 
and finally get the perturbative expansion (L 2 ) = (L (1)2 ) + (L (2)2 ) +2 (L (1) -L (3) ). The term (L (1) -L (2) ) is zero for an underlying 
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Figure 2. Wave vector parts of higher-order corrections to {L 2 } computed for power law spectra P^{p) oc p n , filtered with a Gaussian 
smoothing function at R = 0.5 ft -1 Mpc, as a function of spectral index n, in units of the square of the mass variance <j 2 in order to 
eliminate the dependence on the normalisation of the spectrum. The points corresponding to the CDM spectrum are indicated by stars 
and positioned arbitrarily at n = —1.5 which resembles the slope of CDM power spectrum at galactic scales. 



Gaussian matter distribution, but it should be taken into account in the framework of more general non-Gaussian statistics 
(work in progress). Assuming Gaussian statistics here, we disregard it. In sections 2.1 and 2.2 (for the Einstein-de Sitter 
universe; in Appendix A for a more general Friedmann universe) we reviewed the Lagrangian theory and the perturbative 
solutions Si, S2 and S3 of the Lagrangian fluid equations. Using these results we calculate the corrections (l/ 2 - 12 } and (L^-l/ 3 '} 
in section 3, after summarising the results of linear theory (i.e., the term (l/ 1 -' 2 }). The final expressions are rather cumbersome 
(the details of the calculations have been deferred to Appendix B), but we can summarise the main features of our results as 
follows: for an Einstein-de Sitter universe, 

• (I/ 1 ) 2 ) 1 / 2 ocr" 3 oct [oc D(t)] ■ 
. (L( 2 > 2 ) 1/2 ocr- 5 oc^ 3 [oc-E(r)] ; 
. (L«-Lf ) 1 / 2 ocT- s oct 5 / 3 

• (L^-I/ 12 )) 1 / 2 ocr" 5 oct 5 / 3 

where D is the growth factor of the density perturbations, E and Fh (h = a, b, c) are the growing modes of the second- and 
third-order Lagrangian displacements, respectively. We see that the perturbative corrections to (L' 1 ' 2 } 1 / 2 grow proportionally 
to t 5//3 in the Einstein-de Sitter universe, in agreement with Peebles (1969). The expressions between square brackets give 
the generalisations of the results for a generic Friedmann universe. Interestingly, all the corrections we have analysed are 
proportional to the same invariant of the inertia tensor J of the matter contained in the homogeneous Lagrangian volume V, 
a result we can express as 

. (L^ 2 ) oc <L< 2 > 2 ) oc (L (1) ■ L (3) ) oc /I 2 — 3/i 2 , 



cc (D(r)F h (r)) 



oc (t)(r)[D(r)E(r) - D(r)E(r)]Y 
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Figure 3. Wave vector parts of non-linear corrections to (L 2 ) for a standard CDM spectrum as a function of filtering scale R. Left 
panel: scaling of E< 2 >, E< 3) , E^ 3 ' and E< 12 ) versus R, arbitrarily normalised. The run of the mass rms a with R is shown for comparison 

(same arbitrary normalisation). The values for E^ and Ej. 3 are practically indistinguishable on this plot. Right panel: the various E's 
normalised to cr 4 as a function of filtering scale R. These ratios are independent of the normalisation of the power spectrum. 

where fii and fi2 are the first and the second invariant of the inertia tensor (see equations (40) and (41)). This invariant 
Hi — 3/i2 has been thoroughly investigated in Catelan & Theuns (1996). As a consequence of this factorisation we have been 
able to express the order of magnitude of the non-linear corrections to (L 2 ) in terms of the linear contribution, (L 2 ) = 
(1 + T) (L/ 1 ' 2 } (equation (54)), where T w 0.6 for the standard CDM spectrum at galactic scales. Taking into account that 
the non-linear correction is small, we conclude that linear theory gives a good description of the angular momentum up to 
maximum expansion. Since in addition linear theory predicts, in the Einstein-de Sitter model, a growth rate L oc t, it follows 
that the initial torque is a good estimate for the tidal torque over the whole period during which the object is spun up: 
dL(t)/dt « dL(0)/dt. 

Finally, as is the case with almost any analytic calculation, comparison of these results against observations is hampered 
by the fact that the very final stages of galaxy formation are likely to be highly non-linear and in addition dissipative processes 
may play an important role as well. Analytic investigations are not able to take such highly complex phenomena into account. 
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APPENDIX A 

In this first appendix we report the expressions of the Lagrangian perturbative solutions S n explicitly up to third-order and 
valid for a generic non-flat universe. Specifically, we give the expressions for the growing modes D, E, F a , Ft, and F c for linear, 
second- and third-order terms and the expressions for their wave vector dependence, i.e. the Fourier transforms of the longitudi- 
nal potentials ip^ 2 \ ip!^\ and of the transverse components T a . For brevity we use the symbol 6(r) = ln^r - l)/(r + 1) 
for the open universe case (k — —1) and A(r) = arctang(l/r) = — iQ(ir) for the closed universe case (k — +1). 



A.l Zel'dovich approximation 

The growing mode of the density fluctuations D(t) is given by (Shandarin 1980): 

D{r) = \ |l + 3(r 2 -l)[l + re(r)]} , (60) 
for the open universe, and 

D{r) = \ {-1 + 3 (r 2 + !)[!- tA(t)]} , (61) 

for the closed universe. The solution (61) can be obtained from (60) by substituting in the latter r by it and reversing the 
sign to make the growing mode positive. Note that, in contrast to Bouchet et al. (1992) and Catelan (1995), we normalised 
D(t) according to the suggestion of Shandarin (1980): the coefficient 5/2 is such that D(r) — > t~ 2 in the limit r — > — oo, 
which coincides with the Einstein-de Sitter case. D(t) for the different universes is plotted in Fig. 4. 



A. 2 Second-order approximation 

The time dependence of second-order growing mode E(r) corresponding to the normalisation chosen for D is: 

E{T) = ~f -j(r 2 -l){l+r9(r) + i[r + (r 2 -l)e(r)] 2 } , (62) 

for the open universe and 

£ M = -f + y(r 2 + l){l-rA(r)-i[r-(r 2 + l)A(r)] 2 } , (63) 

for the closed universe. These solutions have been first derived by Bouchet et al. (1992). The extra factor 25/4 of the present 
version is due to the different normalisation of the first-order solution D. An excellent approximation of the second-order 
growing mode is E = — |D 2 (see Fig. 4 for a plot of E(r) and Fig. 5 for a plot of the approximation). In the limit r — > — oo 
one has E = — §t~ 4 which corresponds to the flat case. 
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The Fourier transform of the second-order potential ip( ' is (Catelan 1995): 



1 f d Pi d 



) 6 



(27r) 3 < 5 D (p 1 + p 2 -p) 



K (2) ( Pl , P2 )^( Pl )V il; (p2) 



where we have defined the kernel 



(2) I \ — 1 

K X ''(pi,P2) = ^ 



2 2/ n2 
PlP2 - (Pi - P2) 



(64) 



(65) 



This kernel describes the effects of non-linearity on the second-order Lagrangian motion of the mass fluid elements. 



A. 3 Third-order approximation 

The third-order solution S3 actually corresponds to three separable modes, two longitudinal and one transverse, as was 
discussed previously (section 2.2). The expressions of the growing modes F a ,Ft and F c for a non-flat universe are known up 
to quadratures in terms of the lower-order solutions D and E. Explicitly (Catelan 1995): 



F a (r) = -2D(t) I dnD(ri)- 2 j % dr 2 q(t 2 ) D(t 2 ) 4 , 

J — OQ J —OQ 

F b {r) = -2D(r) f dn Din)- 2 j % dr 2 a(r 2 ) D(r 2 ) 2 [E(t 2 ) - D(t 2 ) 2 

J — OQ J — OQ 

/T f>Tl 
dn / dr 2 a(r 2 ) Z>(r 2 ) 3 . 
- oc J — OQ 



(66) 
(67) 
(68) 



Excellent fits of these functions are F a (r) = — ^D 3 ,Fb(r) = ^f-D 3 and F c (r) — —jD 3 (see again Fig. 4 for a plot of the 
functions and Fig. 5 for a plot of the approximations; these functions were not shown in Catelan 1995). In the limit r — > —00 
one recovers the flat solutions as required: F a (r) — — ^t^ 6 , F b (r) — |yT~ 6 and F c (t) = — i-r -6 . 



The Fourier transforms of the third-order potentials ip^ and rp^ 1 and of the transverse components T a are respectively 
(Catelan 1995) 



/,(3) 



^i 3) (p) 
^ 3) (P) 

T< s >(p) = i 



_ J_ f rfplrfp 
P 2 J (2n 



^[(2^) 3 fo(p 1 +p 2 + p3-p)]^ ) (p 1 , P2 , P 3)^ (1) (pi)^ (1) (P2)^ (1) (p3), 



; /' rf Pl dp 2 dp 3 [(27r)3fo(pi + p2 + p3 _ p)] K (3 )(pij pa> ps) ^(i )(pi) ^(i )(pa) ^(i )(p3)) 



(27r)» 

' /P ^ P :: /1>3 1 + P ,+p: ! -pm;; !; fp, P ,.p: ! )r'-(p 1 .) ( - "(p 2 )r' Jl (p:,). 



where we have introduced the following kernels: 

«4 3) (Pi , P2 , P3 ) = \ e a7 J e/3r,<7 Pa Pl/3 P2 7 P217 P35 P3tr , 



(Pl,P2,P 3 ) 



<4 3) (Pl,P2,P 3 ) = 



P Pi 



Pi - (P2 +P3) P- (P2 +P3, 



|P2+P3| |P2+p 3 | 



K (2) (p 2 ,p 3 ) , 



1 Pi - (P2 + P3) 



K (2) (P2,P3) (P2 +P3 - Pl)a 



(69) 
(70) 
(71) 

(72) 
(73) 
(74) 



2 |P2+P3| 2 

These expressions are more suitable to compute correction terms to L than the original equations as derived in Catelan (1995) 



APPENDIX B 

In this appendix we give explicit derivations for the ensemble averages (l/ 2 ) 2 }, (i/ 1 ) • L^ 3 - 1 } and (L/ 1 '- I/ 12 '}. We will use the 
perturbative corrections to the Lagrangian displacement as reviewed previously. 

• Starting with the former one, we use the expression (26) for L (2) to find: 

<L< 2 > 2 > = E(r) 2 e aM e a9W J al J a , Y { V^Vf^) . (75) 
The second-order deformation tensor 2?^ may be written in terms of the second-order potential ip^ as 
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Figure 4. Growth factors of the first, second and third-order Lagrangian displacements, respectively D(t), E(t) and F a (r), for open 
(upper curves), flat (middle curves) and closed (lower curves) universes. The remaining third-order growth factors Ff>(r) and F c (t) are 
not shown as they are almost identical to F a (r). 

Vt] ee d a d^ {2) (0) = - J ^ Pa Pfj ^ 2) (p) W(pR) , (76) 

where the field is now assumed to be filtered on scales R using the smoothing function Wr. Inserting the expression (64) 
for the Fourier transform i/>' 2 )(p) of tp^ 2 - 1 we find 

f d ^rw (P1 + P2) l (Pl | 2 + P2h + n\ R ) * (2) (p^> p=») ^ (1) (p0 ^ (1) ( P2 ) ■ (77) 

p J (27r) b IPi+PaT 
From this last expression we obtain: 

/<n(2)<n( 2 ) \ / dpi C?P2 dp'i dpa ~,, Id\TI//I ' i '|n\ (2)/ \ (2)/ / / >, 

(^j^v) = / ( 27r )6 ( 27r )6 W(\pi+p 2 \R)W(\p 1 +p 2 \R)K y >(p 1 ,p 2 )K y '(p^p-i) 

(^ (1) (p 1 )^ (1) (p 2 )^ 1) (p' 1 )^ tl, (P2)) ■ <7S! 



(Pi+P2)^(pi + P 2 ) CT (pi + Pa)/s' (Pi + PaV , -.I..,, , -i'W>.V"i«' 



|Pl+P2p |Pi+P 2 | 2 

Since the primordial gravitational potential is assumed to be Gaussian distributed, one has (see, e.g., Peebles 1980) 

<^ (1) ( P i) ^ (1) (p 2 ) Vi (1) (pi) ^ (1) (P2) ) = (2^) 3 fo(px + P2) Pi,(px) (2tv) 3 Snip', + Pi) ^(pi) 

+ (2tt) 3 <J D (pi + pi) P^(pr) (2tt) 3 5 d {P2 + p 2 ) P/>(p 2 ) 

+ (2tt) 3 fo(pi + p' 2 ) P^( Pl ) (2tt) 3 S D (p 2 + pi) P^(p 2 ) . (79) 

The first term does not contribute to the integral since ft^(p, — p) = . The remaining two terms give identical contributions, 
hence: 

,(2W2) , _ » [dpi dp 2 (pi + p 2 )/a (pi + p 2 ) CT (pi + ggV (pi + p 2 )„' r (2 ) , n „^[wnn.xnjmi 2 i 



<^^> = 2/^^ — — 7„ |4 [« W (Pi,P«)] J [^(IPi + P2 |P)]^(Pi)^(P2). (80) 

J \^ ) IP 1 i P 2 

The trick now is to reduce this integral in such a way that we can apply the rule, 
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Figure 5. Ratios between the higher-order growth factors E, F a , Fi, and F c and their approximations in terms of powers of D, for open 
(upper curves), flat (middle curves) and closed (lower curves) universes. The third-order ratios have been offset artificially by 0.1, 0.2 
and 0.3 for clarity. The second-order ratio E/(—3D 2 /7) is shown here for completeness; it has been first shown as a function of fl in 
Bouchet et al. (1992). 



/ dpp a pf3p~,p s F(|p|) = jg{S a p S~,s + S a -, S/3s + S aS Sp-f) / dpp G F{p) , (81) 

J sphere J 

which holds for any function F(p) which depends only on the modulus p of p. Let us define the new variables ki = pi + p2 
and k2 = pi — P2. The determinant of the Jacobian of this transformation is 1/8. At this point, noting that 

„9 1 2 



^P 2 - 4 [^( Pl) P 2 )f=4_ (&i& _ 1 + fc2fc _ 1)2 _^_ 
where 8 = ki - k2/fci k 2 and ku = \kh\, the integral (80) may be written as 



(82) 



dki dk2 kip k\ a kxpi k lo 



(ki k^ 1 + k 2 fcf 1 ) 2 - 40 2 



[W(ki R)] 2 P 4 (|ki + k 2 |/2) P 5 (|ki - k 2 |/2) . (83) 



This expression holds for any power spectrum Ps(p) = p 4 P^(p), but let us assume for simplicity that the power spectrum is 
scale-free, P^>(p) = Ap n ~ 4 . The calculation can be continued as follows. Since one has: 

|ki + k 2 | n |ki - k 2 |™ = (*i k 2 ) n [(fci kz 1 + k 2 fcr 1 ) 2 - 46? 2 ] n/2 , 
the integral in (83) may be simplified to 



(84) 



/t>( 2 )t>( 2 ) \ — : 

\ U 0<j U 0' a' I ~ 



A 2 f dki kip k lrT fc lj3 > fc la / 



(2tt) 3 kf 



[W(ki R)] 



o (2t) 



dk 2 fc „ +2 



2 0/2 



+ 1 



d6 (1- 



fki fey 
U 2 fcJ 



1 (n-4)/2 



-40 1 



We note at this point that the function 

poo /> + l 

2 2„-l A (2) (fci;n) _ J dk2k n+2 J d6 (1-6 2 ) 2 [{k!^ 1 +k 2 K 1 f -Ae 2 }^' 2 
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which typically has to be evaluated numerically, depends only on the modulus |ki| and on the spectral index n, as indicated. 
We can therefore apply the rule in (81) since the smoothing function depends only on the modulus of ki as well, to do the 
transformation: 



dk 



■ hp h a V K' [W(kR)f k n ~ 4 A (2) (pi; n) 
(Spa fip'v' + Sep* 8 aa , + 8p a i 8pi a ) £ (2) (R; n) , 



(87) 



2(2tt) 2 J (2tt) 3 
1 

15 

where S' 2 - 1 is given in equation (43). Finally, following along the lines of the the derivation of (l/ 1 ' 2 } reported in Appendix A 
of Catelan and Theuns (1996), one ends up with the result (42) in the main text. 

The last expressions are not valid for more physical, non-power law, spectra. However, the appropriate expressions can be 
found by following the same strategy. Let us consider for example the case of the CDM power spectrum (see, e.g., Efstathiou 
1989), where the power spectrum is parametrised by: 



P^(p) = Ap~ 3 [T(p)} 2 = Ap- 3 l + (ap+ (bp) 3/2 + (cp) 



-2/v 



where a = 6.4 (Qh 2 )' 1 Mpc, b = 3.0 (ft/i 2 )" 1 Mpc, c = 1.7 (n/i 2 ) -1 Mpc, and v = 1.13; A is the normalisation of the 
spectrum, as before. In this case, one obtains instead of equation (85): 

a2 



4 J (2tt)3 



fa/3 kla fa/3' ki a > 



kf 



[W{ki R)Y 



dk2 



/o (2^) J-i 
where we have defined the transfer functions 



(l - e 2 f [{ki k 2 x + k 2 fa" 1 ) 2 - 40 2 ] 3/2 [T(fa, fa; +9)] 2 [T(fa, fa; -6)] 2 



T(ki,k 2 ;±9) 2 



1 + 



iG{±8) + (bG(±6)y 2 + (cG(±0)) 



i»1 -2/k 



(89) 



(90) 



and G(±8) = ±\/fafa~(fa k^ 1 + k 2 fa" 1 ± 2(9) 1/2 . As a consequence of the change of variables, (pi, P2) — > (ki, k2) we have, for 
example, T(ki,k 2 ;6) = T(pi). Finally, proceeding as earlier we get: 
.2 



<2) 
-'CDM 



(R) 



A 



(2-) 4 J 



dpp*[W(pR)] 2 A (2) DM (p), 



v (2) 
*CDM 



(*l) 



1 f°° f + 1 

-J dkikl J d9(l-0 2 ) [(fa fa 1 + fa fa: 1 



) 2 - 48 2 ] 3/2 [T(fa, fa; +9)} 2 [T(ki, fa; -6)f 



(91) 
(92) 



• Let us now summarise how to deal with the simpler case of the averages (L^ ■ L^}. We first show that (L^ 1 - 1 - Li 3 ') = 0. 
One has: 



<L«- Li 3 > ) = ^(Li 1 )^) = Z)(r)F (r)e a/37 e a/3V J„, J a ,,, (V^V^,J 



(93) 



Writing T> L , as 

2>< 3 > v , ee fy, a CT ,^ 3 > (0) = - J ^ p„, Pu 1 i>i 3) (p) ^(pi?) 
and inserting the expression (64) for ' (p) , one gets 



(94) 



/ dpia 
127 



) 



6 P 2 [W(P2-R)] P^{pi)Pii>{$2)V20'P2o<P2<x"V2fiV2o 



Pl/3" Pl 7 " Plr," P2.5" P2<x" + Pl/3" P2 7 " P 2 r," P«" Pier" + P2/3" Pl 7 " Plr," PlS" Pier" 



(95) 



which is zero for any and power spectrum P^, since symmetric tensors saturate antisymmetric tensors. This completes the 
proof. 



• Next, let us compute: 

(L (D. L (3) ) = ^{L^Lfl) = D(r)F b (r)e a01 e a ^ J„, J^, {V^V^,) 



(96) 
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Non-linear evolution of the tidal angular momentum 17 

The third-order deformation tensor 2?^^ may be written as 

Vfl p ee d a d,^ (p) = -Jjto sPa P0 ^ (p) W(pR) , (97) 
where the field ip^ is now filtered on scale R as previously. Inserting the expression (65) for (p); we get: 
V ba0 = /— (2^9 | Pl + P2 + P3P H/(| Pl +p 2 + p 3 ji?) 



1 2 

X oPl 



1 - 



Pi - (P2 + P3 
Pl|P2 + P3 



K (2) (P2,P3)^ <1) (Pl)V ;(1) (P2)V ;(1) (P3) • (98) 



The average can now be computed easily by following the same steps as in the derivation of (l/ 2 ' 2 }. One ends up with the 
expression 



Wc stress that this results holds for any spectrum P$(p), including the CDM spectrum (88). Limiting ourselves to a scale- free 
power spectrum P^(p) = Ap n and applying the rule (81), one recovers the equation (46) in the main text. The generalisation 
to a CDM power spectrum is straightforward, since there is no coupling between different wave vectors in the power spectrum 
in this case. 

• In a similar fashion, one computes the average 

<L« L< 3 > ) = ^(L«4 3 )) = D(T)F c (T)e a ^e ag , Y J^J^ {V% T%) . (100) 



From the expression (71) for T"i 3 ^(p) we have: 



- 'J dP \Zy PA (Pl +P2 +Ps)f3 ^ 3 V'P2>P3)^(lPl +P2 +P3|^)^ (1) (pi)Vi (1) (p2)^ 1) (p3) ■ (101) 

One now proceeds as earlier. For any power spectrum P^(p), the result is: 

Assuming a scale- free power spectrum P^(p) = Ap n and applying again the rule (81), one recovers the equation (47) in the 
main text. Once again, the generalisation to a CDM power spectrum is straightforward. 

• We conclude this appendix by discussing the third-order correction L/ 12 ^: 



L^^(t) = Vo (DE - ED) J dqV# ' x X7ip [ z) . (103) 
Expanding the potentials in Taylor series around q = 0, one obtains for the a-component 



-(12) 



(DE - ED) e a/37 



Vo F (0) (0) + V™ T>% J*, = + ^ • (104) 



(12) . r(12) 



In general, the first of these two terms may be written as 

-<12) _ ., , :,./:• [ d Pl dp 2 d P3 (P2+P3) 

(2tt)9 Pl/3 |p 2 + p 3 | 

X ^(P2,P3)^ (1) (P1)^ (1) (P2)^ (1) (P3), (105) 



L<™ = noV(DE-Eb)e^ / pl0 ^^W(p,R)W(\p 2 + p 3 \R) 



and the second term as 

r(12) /ni^ z^tS\ f dpidp 2 dp 3 (p 2 + P3) 7 (P2 + p 3 )^ d\w/i im 

l b1 = -(DE-ED)e a01 J a7] J ^-^ JWio- - |p 2 + p 3 p W fa R ) + P*\ R ) 

X K (2) ( P2 ,P3)^ (1) ( P1 )^ (1) (P2)^ (1) (P3) ■ (106) 

With these informations, one can compute the ensemble average (L' 1 ^ • L' 12 ') = (L 1 - 1 ^ • L^" 2 '} + (L' 1 ^ ■ L^ 2 '}. The resulting 
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expression are, respectively, 
. ' ! ■. i. _„..r h,nr- / h) J (i P' f/ P 

|pi - P2|" 2 K (2) (pi,p2)p2c«(pi Q P2/3 -P2aPlf3)p2 1 Jf}~, 

-TjoT T \pia(j>\aP'2[i - Pl<y.P\(3)Pl~, J/3-y] (107) 



(L W -LD = -VorD(DE-ED) \ -^±W(p 1 R)W{p 2 R)W{\Y> 1 - p 2 |J?) P^(pi) P^(p 2 ) 



(L (1) .L^ 12) ) = D(DE-ED) \ ^l^W{p,R)W{ P 2R)W{\p 1 ~ p 2 \R) P^{pi) P^) 



(i).t.( 12 )\ - h/nR.Km f^E* 

(2tt) 6 

X Pi — P2p 2 K (2) (pi,P2)p2a(pi a P2/3 ~ P2a£>l/3)jJ2 7 J"/3 7 (Pli(P2 - Pl)r, Ji,) 

= J r [p2c t (piaP2/3 - P2a£>l/3)P2 7 J/3 7 (Pl«(P2 - Pl)r, Js V )] , (108) 

where we have introduced the integral operator T for conciseness. We proceed by showing that the first term (L/ 1 '- L^ 12 ^} 
is zero for any power spectrum. In the eigenframe of the inertia tensor and taking advantage of the fact that T is a linear 
operator, we find: 

(L (1) • L^ 12) ) = -77 r i a T\A a \ (109) 

where we defined A a = (pi • P2)p| Q — P2PiaP2a- Using the fact that, because of rotational symmetry, f[A ] = F[Ap] for any 
Cartesian axes the result simplifies to 

<L (1) -l4 12) > = - VQ r^T[A z ] . (110) 

Since T, a A a = we finally get JFJEcA^] = = g^A*], which completes the proof. 

We can use a similar trick to simplify the second term: defining B a = pi a (p2 — Pi)a (no sum over a!) it is straightforward 
to show that in the eigenframe of J one gets 

<L (1) -l4 12) ) = r[L a A a i B B ] 

~ ^[L X AxB X "I" LyAyBy + L Z A. Z B Z + lxLy(A X By + A_y B X ) + L X L Z (A_ X B Z + A_ Z B X ) ~T~ LyL Z (AyB Z + .Ag-By)] 

= (/i 2 - 2 M2 )J r [A z B z ] + 2 t i 2 J r [A z B x ] 

= (nl - 3i*2)F[A,B,] . (Ill) 

The last equality is obtained as follows: we first use the same trick as before to prove that f[E B« A z ] = ±F[E a B a XgAp] = 0. 
Next, since because of symmetry J r [A z B x \ — J r [A z B y ] we get J-[A Z B X ] — — }^T\A Z B Z \. The integral T\A Z B Z \ has been 
computed numerically. As in all previously discussed correction terms, \x\ — 3fj,2 factors out. One recovers equation (53). 
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